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We examine a proposal by Sherson and M0lmer to generate polarization-squeezed light in terms of 
the quantum stochastic calculus (QSC). We investigate the statistics of the output field and confirm 
their results using the QSC formalism. In addition, we study the atomic dynamics of the system 
and find that this setup can produce up to 3 dB of atomic spin squeezing. 
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I. INTRODUCTION 

Many recent experiments in quantum information pro- 
tocols and precision metrology have utilized the optical 
Faraday rotation of light passing through spin polarized 
atomic samples. Of particular interest to our discussion 
is a recent article by Sherson and M0lmer, in which it is 
shown that polarization squeezed light can be generated 
by sending a cw beam or pulse of linearly polarized light 
through an atomic gas twice in different directions (See 
Fig. 1) pp. For a proposal involving a similar experi- 
mental setup in the context of atomic magnetometry, see 

0- 

The aim of this paper is to analyze this experiment us- 
ing the quantum stochastic calculus (QSC). After a brief 
introduction to QSC, we describe the system in terms of 
a quantum stochastic differential equation (QSDE). Us- 
ing this model, we derive the atomic and field dynamics 
of the double-pass system, confirming their results, and 
furthermore, showing that the resultant atomic states are 
spin-squeezed up to 3 dB. 



II. THE MODEL 

We consider an atomic gas in interaction with the 
quantized electromagnetic field. (Our development of the 
model closely follows the exposition in [9 , where a single 
qubit, rather than an atomic gas, is studied). As the in- 
teraction between the gas and the field is symmetric, we 
can characterize the atomic ground state by the collective 

spin vector J. Furthermore, assuming that the atoms are 
strongly spin polarized in the x direction, we can apply 
the Holstein-Primakoff approximation and introduce the 
dimensionless variables (x at ,p a t) — {Jy, Jz)/y/Jx- Set- 
ting h = 1, these obey the canonical commutation rela- 
tionship [XauPat] = i- 

The atomic gas is described by the Hilbert space of 
quadratically integrable functions L 2 (R) and the electro- 
magnetic field by the symmetric Fock space T over L 2 (R) 






FIG. 1: Sherson and Maimer's setup for the generation of 
squeezed light. After transmission through a gas of spin po- 
larized atoms, the polarization of the cw beam or pulse of 
linearly polarized light is rotated with a quarter wave plate, 
and it is then transmitted through the gas a second time from 
a different direction. 



(Hilbert space of single photon wave functions), i.e. 



T 



fc=i 
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The Fock space T allows us to characterize superposi- 
tions of field-states with different numbers of photons. 
The joint system of the gas and field together is given by 
the Hilbert space L 2 (R) ® T. 

We examine the interaction between the gas and the 
electromagnetic field in the weak coupling limit [3] , [4] , [5] 
so that in the interaction picture, the unitary dynamics 
of the gas and the field together is given by a quantum 
stochastic differential equation (QSDE) as described by 
Hudson and Parthasarathy [6]. 



dU t = {-L*dA t 



LdA* 



\L*Ldt 



iHdt}U u (1) 
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H = ^a 2 (px + xp), and Uq = I. (For 
a derivation of the double-pass QSDE, see Appendix [A| . 
Here, the operators L and L* are proportional to the 
the annihilation and creation operators on the atomic 
gas system and A t and A^ denote the field annihilation 
and creation processes. Note that the evolution Ut acts 



nontrivially on the combined system L 2 (R) &J 7 . That is, 
L and A t are understood to designate the single system- 
operators L <g> I and I <g> A t respectively. Equation (fiT) 
should be taken as a shorthand for the integral equation 

U t = I- [ L*U T dA T + / LU T dA% -\\ L*LU T dr 

Jo Jo Jo 



i I HU T dr, 
to 

where the integrals on the right-hand side are Hudson- 
Parthasarathy stochastic integrals, and where the 
stochastic increments satisfy the following rule jB], [7]: 

Quantum Ito rule Let X t and Y t be stochastic inte- 
grals of the form 



dX t = C t dA t + D t dA* 



E f dt 



dY t = F t dA t + G t dA\ + H t dt 

for stochastically integrable processes C tl D tl E tl F tl G t 
and H t (see [6], [7] for definitions). Then the process 
X t Y t is itself a stochastic integral and satisfies 

d(X t Y t ) = X t dY t + (dX t )Y t + dX t dY u 

where dX t dY t should be evaluated according to the quan- 
tum Ito table: 





dA t dA* dt 


dA t 

dA* 
dt 


dt 





i.e. dX t dY t = C t G t dt. 

We now describe a simple notation (introduced in [9]) 
for expressing differentials of products of stochastic inte- 
grals: Let {^i}i=i,...,p be stochastic integrals. Then we 
write 



d(Z 1 Z 2 ...Z p ) 



£ M> 



i/C{l,.-,p} 



where the sum is taken over non-empty subsets of 
{l,...,p}. For v — {ii, . . . ,ifc}, the term {v} is ob- 
tained by differentiating only the terms with indices in 
the set {ii, . . . , ifc} and preserving the order of the fac- 
tors in the product. For instance, consider the differen- 
tial d(Z\Z 2 Z 3 \ which contains terms of type {1}, {2}, 
{3}, {12}, {13}, {23} and {123}. Applying the Ito 
rule, we have {2} = Z 1 (dZ 2 )Z 3 , {13} = (dZ 1 )Z 2 (dZ 3 ), 
{123} = (dZ 1 )(dZ 2 ){dZ 3 ), and so on. (To avoid confu- 
sion with the notation for references, we have used curly 
braces {} instead of the square braces [ ] used in [9]). 



III. DERIVATION OF THE JOINT 
CHARACTERISTIC FUNCTIONS 

In this section, we derive the joint characteristic func- 
tions for the combined atom-field system. Although we 



will not be using the input-output formalism [8], 
present it here for comparison with pQ: 



we 
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dA t + cL4* 
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nut _ TT*„in TT j^out _ d [ U t( A t + A$)U t ] 

X at — U t X at U t dX vh — — 



Jz/yJx dp 



in 
ph 



vZ l = utfSUt 



dp°X 



dA t - dA* t 

d[U*(A t -A*)U t ] 
iV2 



For the sake of simplicity in later computations, we have 
defined x°^ and p°Jf in such a way that [^/f ,Pp/f ] = it. 
As will be apparent later, once we have calculated the 
variances for these operators, we rescale by a factor of 1/t 
to obtain the actual variances for the modes that we are 
interested in (i.e. the squeezed and anti-squeezed modes). 
Applying the Ito rule, and evaluating the products of 
stochastic integrals using the Ito table, we can express 
the above input-output relations in the language of pQ: 



Pph 

dx°f(t) 



P° P h(t) = p%(t) - ax° a ?(t) 



(2) 
(3) 



dt 

dp%\t) 

dt 



*P l p l(t) (4) 

-a(x%(t)+ap°?(t)) (5) 



As an example, we derive the expression for x^(t): 

d[U t *(A t + A* t )U t ] 

X ph ai — r- 

= {1} + {2} + {3} + {12} + {23} + {13} + {123} 

A simple calculation shows that the terms {1}, {3}, and 
{13} sum to zero, and from the Ito table, we see that 
third powers of increments (i.e. {123}) vanish, leaving 
us with {2}, {12} and {23} to calculate: 



{2} 



dA t + dA* 
T ,dA t + dA* t] 



LI 



{12} = du; """ ' ""* u t = u;^=u t dt 

{23} = {12}* = U t *^=U t dt 



Summing the preceding terms, and substituting for x p \, 
we obtain: 

{ 2 } + {1 2 ) + ( 23) = d ^ +f , 1 -<^±a^ 
= x™ h dt + u;ap%U t dt 

=>xff(t) = 4i(t)+a P :r(t) 



For comparison with the results in [1 , a should be 
equated with ft, and as we have ignored damping, r 
should be set to 1. 

Define F[t,k,l] and G[t,k,l] as follows (see, for in- 
stance, [TT]): 

F[t,k,l] :=H$ | U;(J lp ®e iki ^ il )Ut \ v <g> $) (6) 



G[t,k,l] :=(v®$ | /7 t *(e 



"V^)[/ t | v ®$) (7) 



In this notation, F denotes the joint characteristic func- 
tion for * /- * and p, while G denotes the joint charac- 



V2 



-A* 



teristic function for -. / - t and x. While in general, we 
would need to calculate the joint characteristic function 
for all 4 variables, in the particular system we are study- 
ing, the function decouples into two independent compo- 
nents. Here, F and G are expectation values taken with 
respect to an x-polarized spin state of the atoms and the 
vacuum state of the field, as described in the previous sec- 
tion. Since we are interested in obtaining joint character- 
istic functions, and not individual moments, the atomic 
and field operators are given by complex exponentials of 
the respective observables. 

We can calculate F[t,k,l] and G[t,k,l] by solving par- 
tial differential equations given by the following lemma: 

Lemma: 



~dt 
dG 



= -\(al-k) 2 F-a(al-k)^ 



= -\(al + k) 2 G- 



ak^ 



(8) 
(9) 



where F[0,M] = G[0, Jfc, I] = e~T. 

Proof: Let F(Z) := (v <8> $ | U*(Z <g> e ik ^vs* )U t \ 
v <g) $) so that F[e ilp ] = F[t,k,l]. Using the notation 
introduced previously, we have the expression dF[e llp ] = 
(v<8>$ J {l}+{2}+{3}+{12}+{23}+{13}+{123} | v®$). 
Applying the ltd rule and noting that the third powers 
of increments vanish, we are left with the following dif- 
ferentials to calculate: 



{i} 

{2} 
{3} 



dU?(e Up ® 



C/ t *(e ap ( 

U*(e ilp ( 



<e k V2 )u t 

ik(dA t +dAl) 



V2 

U V2 



., (A t +A*) 

\edty k —^)u t 



{12} = dU*(e llp ®( 

{23} 

{13} 



u;(< 



Up, 



)dU t 



ik(dA t +dA*) 
ik(dA t +dA*) 



dU*(e 



Up 



V2 

(A t +A*) 
V2 



-\kHty k ^^)u t 

\k 2 dt)e lk V2 )dU t 



ik 



)dU t 



In the above expressions, dU t and dU% are given by equa- 
tion M. The terms dA t and dA\ vanish with respect to 



the vacuum expectation (see reference [7]) and we find 
that {1}, {3}, and {13} give the following: 

(v®$ | {l} + {3} + {13} |v®$) 
= -±F[L*Le Up + iHe Up ]dt 

- lF[e ilp L*L + ie ilp H]dt + F[L*e ilp L]dt 
= F[C(e ilp ))dt, 

where C(Z) is the Lindblad operator given by 

C(Z) = -§{L*L, Z} + z[#, Z] + L*ZL 

Recalling that L = ^^, L* = ^M, and ff = 

\a 2 (px + xp), we can expand the Lindblad term as fol- 
lows: 

F[C(e ilp )]dt = F[-\{L*L,e ilp } + i[H,e ilp ]+Le ilp L*]dt 
= -\a 2 l 2 F[e ilp ]dt - ia 2 lF[pe ilp ]dt 



a 2 l 2 F[t,l,k]dt - a 2 l^^dt 



1^,272 
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In the last step, we used the equality F[pe dp ] = 
—i — qY - . Summing the remaining terms, we have: 

(v<8>$ I {2} + {12} + {23} |v®$) 

= -±fc 2 F[e iZp ]dt - ^F[L*e ilp + e iZp L]dt 

= -|ife 2 F[e^]dt - i|F[^±^e^ + e^^^]dt 



k 2 F[e ilp ]dt + ^F[£, jfc, Z]dt + a* 



.U2 

4 



3F[t,Z,fc] 
dl 



dt 



Collecting like terms, we arrive at the expression stated in 
the lemma. The initial condition is obtained by noting 
that the atoms begin in a harmonic oscillator ground 
state with a 2 = a 2 = \. As the Fourier transform of 

•^at Pat Z 

a Gaussian remains Gaussian, we have that F[0, /c, I] = 

e~^ The derivation for G[t, k, I] proceeds analogously. D 
We then arrive at the following solutions for the joint 
characteristic functions: 



Pat Pat, x p h x ph 



F[t,k,l] = e 

G[t,k,l] = e -+-/ +2 <^ fc;+ <* fe2 ] 



where 



Pat i^ph 

crl. 



= Jd- 

e~ 



2a 2 t\ 



-2a 2 t 



and 



4a 



*x at 

T 2 

x at ,p p h 

2 



l T ( 3J -»-2a a * 



-2e~ ta ) 

4e-° 2 *) 



(10) 
(11) 

(12) 
(13) 
(14) 

(15) 
(16) 

The expressions <j 2 denote the variances and covariances 
of the respective quantities. Recalling that [XpffiPph] = 



Ul + aH) 



-Wt 2 



(y^ 



*■ 



la 4 t 3 



it, we define the normalized modes x p h = x p h/Vt,p p h — 
Pph/Vi, so that [x p h,Pph] — i- Inserting explicit t- 
dependence, we then obtain the following variances for 
the normalized modes: 



a 2 



a p p h 



(t) 

(t) 



\ a 2 

1 



(*) = eM3 + - 



t < h (t) 



\aH z 



4e" a *)(18) 
(19) 



Comparing equations (18) and (19) to the ground state 



variances o 2 ~ (0) = cr~"T0) = i, we see that the double- 

•J^ph v ' Pph v ' £ ' 

pass setup generates an arbitrary degree of polarization 
squeezing, and from equations ( [l2| ) and (15), at most 
3 dB of atomic spin squeezing. Note, however, that 
resultant states are not minimum uncertainty. Although 
it is to be expected that for a linear system an initial 
Gaussian state remains Gaussian, the derivation of the 
joint characteristic function makes this fact explicit. 
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APPENDIX A: DERIVATION OF THE 
DOUBLE-PASS QSDE 

We begin with the familiar single-pass QSDE's given 
by the Faraday interaction Hamiltonian (see for ex. [TO]). 



■■-{■ 
> - {■ 



-tap 



-tax- 



dA t - dA* 

iV2 
dA* + dA t 



WpHt\ul (Al) 



V2 



\a 2 x 2 dt 



U? (A2) 



To derive the QSDE for the double-pass system, we 
first divide up the light into pulses at discrete time inter- 
vals and have each "piece" make multiple passes. Con- 
sider the following: 



IV. DISCUSSION 

We have used a quantum stochastic model to investi- 
gate the dynamics of a system in which a laser beam is 
sent twice through an atomic gas in different directions. 
Using this model, we reproduce Sherson and Maimer's re- 
sults demonstrating polarization squeezing of the output 
field, and in addition, we show that this setup generates 
a maximum of 3 dB of atomic spin-squeezing. The pri- 
mary difference in using the quantum stochastic calculus 
is that the derivation operates entirely in the time do- 
main. While the time and frequency domain methods 
are equivalent in the simple linear systems that we ex- 
amine in this paper, the advantage of quantum stochastic 
models is that they can be used to study non-linear sys- 
tems as well. 



U< 



t+dt 



ui 



dUl 



M\Ml_ txt ..MlTJi 



M\M\_^...Mi 



where i = 1,2 
and M} = {I - 



M} 



{I 

dA^+dAt 



iap( 



dA*-dA t 



) - Wdt}, 



V2 ) 4 L 

t — \-L — iOix{ ai±t ^^ t ) — \a 2 x 2 dt}. The last two 
equalities follow by recursively applying the first and the 
fact that Uq = Uq = I. This expression is equivalent to 
the series product in [12] and [13] . We can now write the 
QSDE for the combined system as follows: 



n 21 
u t+dt 



dUi 



21 



(M?M})(M*_ M M}_ M )...{MgMl) 

(MfMDU? 1 
(MfMDU? 1 - U\ 



21 



Expanding the above product using the ltd table and 
dropping the superscripts, we arrive at the double-pass 
QSDE: 



ACKNOWLEDGMENTS 



dU t = {-L*dA t + LdA\ - \L*Ldt - iHdt}U u (A3) 



The authors would like to thank Luc Bouten for in- 
sightful discussions. G.S. is supported by the ARO under 



with L = ^iP^l^ H 
before. 



^a 2 (px + xp), and Uq = I as 



M0lmer, Phys. Rev. Lett. 97, 



[1] J. F. Sherson and K. 

143602 (2006). 
[2] H. L. Partner and B. D. Black and J.M. Geremia, 

arXiv : quant -ph/ 0708 . 2730 (2007) 
[3] J. Gough, [arXiv : quant -ph/041 106 4 (2004) 
[4] J. Gough, Commun. Math. Phys. 254, 489 (2005). 
[5] L. Accardi and A. Frigerio and Y.G. Lu, Commun. Math. 

Phys. 131, 537 (2005). 
[6] R. L. Hudson and K. R. Parthasarathy, Commun. Math. 

Phys. 93, 301 (1984). 
[7] K.R. Parthasarathy, An Introduction to Quantum 



Stochastic Calculus, (Birkhauser, Basel, 1992). 
[8] C. W. Gardiner and M. J. Collett, Phys. Rev. A. 31, 

3761 (1985). 
[9] Bas Janssens and Luc Bouten, Journ. Phys. A. 39, 2773 

(2006). 
[10] L. Bouten and J. Stockton and G. Sarma and H. 

Mabuchi, Phys. Rev. A. 75, 052111 (2007). 
[11] A. Barchielli and G. Lupieri, J. Math. Phys. 26, 2222 

(1985) 
[12] M. Yanagisawa and H. Kimura, IEEE. Trans. Aut. Con. 

48, 2107 (2003). 



[13] J. Gough and M. James, arXiv:quant-ph/0707.0048vl 
(2007) 



